Effect of long range connections on an infinite randomness fixed point associated with 
the quantum phase transitions in a transverse Ising model 
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We study the effect of long-range connections on the infinite-randomness fixed point associated 
with the quantum phase transitions in a transverse Ising model (TIM). The TIM resides on a 
long-range connected lattice where any two sites at a distance r are connected with a non-random 
ferromagnetic bond with a probability that falls algebraically with the distance between the sites as 
l/r d+a . The interplay of the fluctuations due to dilutions together with the quantum fluctuations 
due to the transverse field leads to an interesting critical behaviour. The exponents at the critical 
fixed point (which is an infinite randomness fixed point (IRFP)) are related to the classical "long- 
range" percolation exponents. The most interesting observation is that the gap exponent ip is exactly 
obtained for all values of a and d. Exponents depend on the range parameter a and show a crossover 
to short-range values when a > 2 — t]sr where i]sr is the anomalous dimension for the conventional 
percolation problem. Long-range connections are also found to tune the strength of the Griffiths 
phase. 

PACS numbers: 75.10.Nr, 05.50.+q, 75.10.Jm 
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The presence of quenched randomness drastically mod- 
ifies the nature of the zero-temperature (quantum) phase 
transitions^^. In contrary to the pure systems, many 
low-energy and low frequency properties of systems with 
quenched randomness are dictated by locally ordered 
"rare regions"—. These active rare regions produce 
stronger effects on zero-temperature transitions than on 
finite-temperature classical transitions^. 

Some of the novel features associated with low- 
dimensional random quantum transitions include acti- 
vated (quantum) dynamical scaling at the quantum crit- 
ical point and the existence of Griffiths-McCoy (GM) 
singular regions 6 -^ where the response function diverges 
even away from the critical point. The existence of the 
above mentioned features well-established in the case of 
one dimensional random quantum Ising models by using 
a novel renormalisation group (RG) technique 5 - which is 
exact in the asymptotic limiti. For d = 1, the RG flow on 
the the critical manifold for strong randomness is towards 
an infinite randomness fixed point where the quenched 
randomness effectively grows stronger and stronger as 
the system is coarse-grained^. Numerical studies using 
an extension of the above RG scheme^ as well as quan- 
tum Monte Carlo studies^ predict a similar scenario of 
IRFP even for random quantum Ising transitions for spa- 
tial dimension d = 2. Also at an IRFP, stronger cou- 
plings always dominate and hence frustration turns out to 
be irrelevant. Therefore, the critical behaviour of quan- 
tum Ising spin-glass and random ferromagnetic models 
are expected to be governed by the same fixed point . 
Long-range spatial correlations of disorder are found to 
enhance the off-critical singularities^. 

Situation becomes very interesting when the quenched 
randomness arises due to dilution or vacancies in a quan- 
tum magnetic mode l 12 d 3 i 14 i 15 i 16 . The problem is imme- 
diately connected to the well-studied problem of the per- 
colation transition of the underlying lattice^. The study 
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FIG. 1: The phase diagram of a dilute TIM in the T — p 
plane. Senthil and Sachdev 1 studied the quantum transition 
along the vertical line which terminates at the multi-critical 
point (p c ,r o ). In the long-range connected case, the phase 
diagram is qualitatively similar. However, p c depends on a. 



of quantum phase transitions in a TIM at the percola- 
tion threshold of a dilute lattice shows the existence of an 
infinite randomness critical point scenario in dimensions 
greater than unityi^. 

The experimentally studied systems LiHoF4 which are 
suitably modeled by transverse Ising spin glass mod- 
els however are complicated because of the presence of 
long-range dipolar interactions^. Moreover in a metal- 
lic system, the presence of long-range RKKY interaction 
is expected to modify the critical behaviour—. Question 
therefore remains whether a non-trivial long-range inter- 
action modifies the exponents associated with an IRFP. 
To address this particular issue, a quantum Ising spin- 
glass model with long-range random interactions, where 
interaction between any two spins separated by a dis- 
tance r.ij is Jij/rfj with random J^ 's, was introduced 2 ^. 
The perturbative RG calculations 21 of the above spin- 
glass model around the upper critical dimension failed to 
locate any stable weak coupling fixed point in the non- 
mean-field region and a RG flow towards the strong cou- 
pling was observed. This immediately points to the ex- 
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istence of an IRFP for the above spin-glass^ . However, 
no quantitative information can be derived regarding the 
nature IRFP in the above spin-glass model. 

In this communication, we study the quantum tran- 
sitions in a novel TIM with ferromagnetic interactions 
defined on a long-range bond dilute lattice. It should 
be emphasized at the outset that in this model there is 
no long-range random interaction, nevertheless, there are 
long-range connection probabilities which effectively sim- 
ulate a long-range interaction that non-trivially affects 
the exponents at IRFP. The Hamiltonian of the model is 

H = -Y,JijS!S*-rY,s? (i) 

where S^s are Pauli spin matrices and T is a non- 
commuting transverse field. The bond Jy is chosen from 
the binary distribution 

P ( J v) = -T^ 5 i J ij - J) + ( 1 - ) 5 i J n) ( 2 ) 

where p (0 < p < 1) denotes the connection probabil- 
ity between the nearest neighbour sites, d is the spatial 
dimension and the parameter a determines the range of 
connection. The parameter a is restricted to non-zero 
positive values to ensure extensivity. The Hamiltonian 
refers to a pure TIM which resides on a lattice where 
any two sites separated by a distance r is connected with 
a ferromagnetic bond of strength J with a probability 
falling algebraically with the distance between the sites 
as p/rfj~ a . Clearly, for smaller values of a long-range 
connections are dominant. 

An infinite randomness critical fixed point is charac- 
terised by three exponents^: At the critical point, the 
smallest energy gap AE is related to the linear dimension 
L of the rare region as — In AE ~ which is a signature 
of activated quantum dynamical scaling. Analogously, 
the magnetic moment of the cluster grows as L^. The 
typical size of such a cluster (the correlation length) di- 
verges as £ ~ 8~ u where 5 measures the deviation from 
the quantum critical point. All the bulk exponents are 
expressed in terms of the exponents <j>, ip and Senthil 
and Sachdev 13 explored the phase transition in the same 
TIM described by the Hamiltonian (1) on a conventional 
dilute lattice where only the nearest neighbour sites are 
connected with a ferromagnetic bond of strength J with 
probability p so that 

P(J t3 ) = P 5{J l3 - J) + (1 - p)5{J lj ), (3) 

where (i,j')'s are nearest neighbour sites. The zero- 
temperature phase diagram of the model in the V — p 
plane is shown in the Fig. 1—. For T = 0, the system un- 
dergoes a geometrical phase transition at the percolation 
threshold p c . For p > p c , at a sufficiently large Y = T c (j>), 
the quantum fluctuations destroy the long-range order. 
However, even at p c there exists an infinitely connected 



cluster with fractal dimensions rf/(< d). Even for a one- 
dimensional TIM, the critical transverse field is finite and 
thus a finite amount of transverse field is required to de- 
stroy the long-range order at p — p c . This leads to the 
existence of a multi-critical point with T = r o and p = p c 
as pointed out by Harris 12 . Senthil and Sachdev studied 
the transition along the vertical line (p = p c and T < T a ) 
where the quantum transition is dictated by the percola- 
tion fixed point p c . The quantum critical exponents are 
solely determined by the classical percolation exponents. 

To comprehend the nature of IRFP associated with 
the quantum transition of the TIM defined through 
Eqs. (1) and (2), it is necessary to explore the associ- 
ated percolation transition with long-range connection 
probabilities 22 . To study the critical behaviour of the 
above percolation model, the mapping of the same to the 
corresponding g-state Potts Model (in the limit q — > I) 23 
turns out to be useful. The resulting Hamiltonian is 

ij 1 ij \a = l / 

where S"' s are the g-state Potts model variables and the 
interaction if is a function of nearest neighbour connec- 
tion probability p. Here, the connection probability be- 
tween two sites gets translated into the interaction be- 
tween the spins as we have only retained most dominant 
(or relevant) long-range interacting term in (4). 

Long-range interacting ferromagnetic Potts Model de- 
scribed by Hamiltonian (4), exhibits a non-trivial order- 
disorder transition at a finite temperature T c {a) (even 
when the spatial dimensionality d = 1 if the range pa- 
rameter a < l 22 i 24 ) . The marginal case d = a = 1 shows 
a topological transition as in the inverse-square Ising 
mode l 16 ! 25 . In the corresponding percolation model (with 
connection probability p/rf^) therefore, an infinite per- 
colation cluster always exists if the nearest neighbour 
connection probability p is greater than p c (a) and there 
is a percolation transition at p = p c (o~). For p < p c (cr), 
infinitely connected cluster disappears. The percolation 
threshold p c (cr) shifts to the lower value of p as the range 
parameter a decreases, i.e., the connection probability 
between two distant sites gets lower. At p c {(j), how- 
ever, there exists an infinite cluster with fractal dimen- 
sion df(d > df > 1). Using the same argument as given 
above, we conclude that the phase diagram of the TIM 
on a long-range connected lattice is qualitatively similar 
to that of the conventional dilute magnet^ (Fig. 1). A 
similar phase diagram is also expected for the long-range 
connected TIM even for d = 1 if a < 1—. 

In the same spirit as in reference [13], we shall now ex- 
amine the nature of transition at the percolation thresh- 
old p c (c) with r < r o (Fig. 1). Along this vertical 
line, the disorder averaged static correlation function be- 
tween two spins separated by a distance x ( which is 
actually the probability that two sites separated by a 
distance x belongs to the same cluster) is given as^£ 



3 



C{x) ~ i/( x ( d - 2 +VLR{<r))^ w here r] LR {a) is the anoma- 
lous dimension exponent of the percolation transition 
with long-range connections. This algebraic fall shows 
that the vertical line is critical with transitions dictated 
by the percolation fixed point p c {a). 

To derive the percolation exponents, let us use the 
Potts model Hamiltonian (4) in the q — ► 1 limit (where 
the order disorder transition in the Potts model is always 
continuous). In the continuum limit, the transition in the 
above Potts model is described in terms of a cubic model 
with a long-range interaction term. The corresponding 
Landau-Ginzburg- Wilson action is2£ 



i [■ i 

U = -- (r + ak 2 + bk a )Y,Q tl (k)Qu(k) 
1 J i=i 

>J J J^2Q ll {ki)Q ll {k 2 )Q tl {k 3 )5{k 1 +k 2 + k^) 
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first order in e(= 3a — d), the exponent v p = l/a + e/Aa 2 . 
The exponents not only depend on the range parameter 
a but also cross over to the short-range value as a ap- 
proaches 2 — r/sR when the long-range interaction term 
bk" becomes irrelevant. In the following, we shall use the 
standard percolation scaling relation for the probability 
that a site falls in a finite cluster of size s given as 17 



P( S ,p) = S 1 - T f( S (p-p c (a)) 1 /n 



(6) 



m(h) 



where Qu is a diagonal traceless q x q tensor and J de- 
notes the integration over k in the first Brillioun Zone, r 
and ii are the mass and the coupling terms respectively. 
The term bk a arises due to the long-range interaction 
in the Hamiltonian (4). At the Gaussian level this term 
dominates over the short-range term ak 2 as long as a < 2. 
For higher values of a ( > 2), long-range interactions is 
irrelevant and therefore exponents of the Potts model are 
of short-range nature. In higher-orders when the effect 
of the coupling term w is perturbatively included, the 
short-range term ak 2 picks up the anomalous dimension 
VSR- We therefore conclude that the long-range interac- 
tion is relevant as long as a < 2 — t/sr- The most inter- 
esting feature of any long-range interacting spin model, 
described by action, (5) is that the term bk a does not 
get renormalised at any order of perturbation and hence 
the anomalous dimension for the long-range interacting 
system sticks to the value t)lr = 2 — a for all values of 
a < 2 — t]sr—. It can also be checked by simple dimen- m{h) 
sional counting of different parameters of action (5) that 
the upper-critical dimension and the range follow the re- 
lation d u = 3cr„, i.e., for a given d if the range parameter 
< c u (= d/3), the interactions are long-ranged enough 
to render mean field theory to be exact. For a = 2 (short- 
range case), the upper critical dimension is as usual 6. 

Let us now restate the above results for the equiva- 
lent long-range percolation problem. As expected, the 
exponents depend on a. In particular, the anomalous 
dimension exponent, r) p (a), of a percolation transition 
is related to the fractal dimension of the lattice at p c 
through the scaling relation^ 7 . df(a) = (d + 2 — r) p (a)). 
In the long-range model, rj p (a) — tilr(<j) = 2 — a which 
immediately leads to the result df(a) — (d + <r)/2. The 
fractal dimension at the percolation threshold is therefore 
exactly known for long-range percolation problem for all 
values of d and a. The other exponents of the perco- 
lation transition can be obtained by perturbative calcu- 
lation around the upper critical dimension™ e.g., to the 



where r = (d + df)/df and a = udf. 

We are now in a position to obtain the critical expo- 
nents related to the quantum transitions at p = p c and 
r < Tq. Let us first look at the exponent ijj. The en- 
ergy gap AE at the percolation threshold of a cluster 
of size L scales as AE ~ exp(— cL d f) where c is a non- 
universal constant. The logarithmic scaling of the cor- 
relation length with the frequency can be easily verified 
also by calculating the dynamical response function us- 
ing the percolation scaling relation (6)^. The exponent 
■0 is therefore given as ip = df = (d + cr)/2. The impor- 
tant features which constitute the essential theme of this 
communication need to be highlighted here: i) the expo- 
nent 4>( <J ) is exactly known for all values of a and d while 
in the conventional case the exponent is exactly known 
only in the spatial dimension 17 , d = 2 ii) ifi depends on 
the range of interaction a and thus gets modified as the 
range parameter a is tuned and iii) ip(cr) undergoes a 
crossover to the short-range value when a > 2 — t/sr. 
To evaluate the exponent cp, let us now use the scaling 
relation of magnetisation in a weak magnetic field h at 
the quantum critical point given asA 



"•<T> 



<f>-d/i/> 



(7) 



where h is a non-universal constant. In the present TIM, 
the magnetisation at p = p c scales as 



2-(d+df)/df 



(8) 



Comparing Eqs. (7) and (8) and using ijj = df and the 
percolation exponents^ 7 - fl p = (t — 2) /a, j p — (3 — r)/er, 
we have cf> — (d — (3 v /v p )/df = 1. 

As in the conventional dilute magnel^, in the present 
long-range case also both the sides of the transition point 
p c are found to be flanked with Griffiths phase with con- 
tinually varying exponents. For example, in the disor- 
dered phase (p < p c (a)), the disorder averaged imaginary 
part of local dynamical susceptibility goes as 



d/z 



(9) 



where in arriving at Eq. (9), we have used the same 
functional form of the percolation scaling function f(x) 
(Eq. 6) as the conventional case^ 7 - with exponents de- 
pending on a. Eq. (9) shows that the paramagnetic 
phase is gapless with a power-law density of states at low 
energy. This power-law singularity leads to off-critical 
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GM singularities, e.g., the average local susceptibility di- 
verges even in the paramagnetic phase as T d ' z ~ 1 when 
the temperature T — > 4 . The dynamical exponent 
z is related to the percolation correlation length £ p as 

z = = (p c — p)^ l/d f . The dynamical exponent varies 
continuously with (p c — p) and diverges at the percolation 
threshold. We immediately conclude that the dynami- 
cal exponent z increases with increasing a leading to a 
stronger divergence of z at p c , indicating an enhancement 
in Griffiths-McCoy singularities. Therefore, the strength 
of the Griffiths phase also varies as a is tuned. 

In conclusion, the transitions in a TIM at the perco- 
lation threshold of a dilute lattice with long-range con- 
nection probabilities provide an ideal situation where the 
effect of long-range connections on the IRFP can be ex- 
tensively studied. The exponents ip and v of IRFP are 
found to depend non-trivially on the range parameter a 
and cross over to the short-range values when the range 
parameter a > 2 — i]sr. However, for lower values of 
c(< fu), when the long-range connections are prominent, 
the exponents are described in terms of mean field per- 
colation exponents 1 "^. Most importantly, the exponent ip 
is exactly known for all values of d and a and increases 
as (7 increases. The strength of GM singularities also get 
enhanced for higher values of a. This is expected because 
for smaller values of a(< a u ) the model becomes more 
mean field like and GM singularities should be weak. A 
recent numerical study (based on the finite size scaling 



argument) of a percolation transition in a power-law di- 
lute chain? 8 in d — 1 suggests that the scaling relation 
df = (d + a)/2, obtained from renormalisation group ar- 
guments, may not be valid when long-range connection 
probabilities are dominant, i.e., for < a < 1/2. How- 
ever, for higher values of a the numerical results corrob- 
orate the renormalisation group results. No such numer- 
ical study is available for the percolation transition with 
power-law dilution in dimensionality greater than unity. 
Question therefore remains whether one should use the 
value of df as obtained from the Potts model analogy in 
smaller a region for d = 1. Also, Vojta and Schmalian 14 
showed that for dilute quantum rotors the quantum dy- 
namics is of conventional (power-law) nature (for tran- 
sitions below the multicritcal point) with the dynamical 
exponent z — df and the correlation length exponent 
v = v p . The scaling relation employed in Ref. [14] 
should in principle also work in the case of the present 
dilute lattice with long-range connections. We should 
also mention that the present work may also be relevant 
for the studies of quantum phase transitions on some in- 
teresting network models 29 . 
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